Integrated photonic qubit quantum computing on a superconducting chip 
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We study a quantum computing system using microwave photons in transmission line resonators 
on a superconducting chip as qubits. We show that all control necessary for quantum computing 
can be implemented by coupling to Josephson devices on the same chip. We take advantage of the 
strong nonlinearities inherent in Josephson junctions to realize qubit interactions. We analyze the 
gate error rate to demonstrate that our scheme is realistic even for Josephson devices with limited 
decoherence times. A conceptually innovative solution based on existing technologies, our scheme 
provides an integrated and scalable approach to the next key milestone for photonic qubit quantum 
computing. 
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Despite the vast potential of quantum computers, no 
perfect physical implementation has been found for quan- 
tum computers. This can be seen by examining two rep- 
resentative systems. Josephson device based supercon- 
ducting systems are easily integrable and scalable, but 
are plagued by the short decoherence times of Josephson 
qubits due to coupling to their complex solid-state envi- 
ronment. Photonic qubits, which have superb coherence 
properties, suffer from the fact that photons do not in- 
teract easily. Also, systems based on conventional bulk 
optical devices are hard to miniaturize and scale. 

Recognizing the importance of integrated systems for 
scalable quantum computing, a number of investigators 
have demonstrated on-chip waveguide based quantum 
gates for photonic qubits recently This is a signifi- 
cant step that may represent the future direction of pho- 
tonic qubit quantum computing technologies. However, 
it is a daunting task to achieve a fully integrated pho- 
tonic qubit quantum computer using conventional tech- 
nologies including those developed in the latest experi- 
ments. This is because conventional optical devices such 
as lasers, lenses, optical cavities and photo detectors are 
bulk devices based on very different technologies and no 
process exists yet to integrate them on the same chip 2] . 
Therefore, alternative realistic approaches to fully inte- 
grated photonic qubit quantum computing systems are 
highly valuable. 

We combine the strengths of photonic and supercon- 
ducting systems to realize fully integrated photonic qubit 
quantum computing. Our physical system is a super- 
conducting chip on which high-Q transmission line res- 
onators (TLRs) and Josephson devices are fabricated. 
The same system has been used for study of cavity QED 
based on Josephson qubits 0, 0, @|. However, in our 
scheme the quantum information is carried by the mi- 
crowave photon modes in the TLRs and the Josephson 
junctions play the role of optical devices. For high-Q 



are available because Josephson devices can be fabricated 
with great precision and controlled conveniently by mon- 
itoring their electrical signals. A further key advantage is 
we can use the strong nonlinearities inherent in Joseph- 
son devices to induce interactions between photons. It is 
shown that high gate fidelities can be achieved even for 
Josephson devices with limited decoherence times making 
their unavoidable noisy environment no longer a limiting 
factor. Therefore, our scheme is a realistic approach to 
scalable photonic qubit quantum computing. 

We start by considering the two identical TLRs shown 
in Figllja) . The TLR mode frequencies are given hy ui — 
n7r/V LC, n an integer and L and C the total inductance 
and capacitance of the TLR. We use the n = 2 mode. For 
L = 0.5nH and C — 5pF, its frequency wq/Ztt « 20GHz. 
The second-quantized voltage and current associated 
with this mode isV{x,t) = y/ hioo / C cos {a{t) + d^t)) 
and I{x, t) — ~ i y^hiuo/L sin ^^{a{t) — a^{t)), where I the 
length of the TLR, x e [—1/2, 1/2] the position along the 
TLR, and d{t) = ae"*'^"* the mode's annihilation opera- 
tor. 
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TLRs the photons have a long life time [I4I which is a 
major advantage. Easy operation and accurate control 



FIG. 1: (a) A photonic qubit based on two TLRs. A and B 
are capacitively coupled to the coupling CBJJ via Cc- B is 
also coupled to the right CBJJ via Cr- (b) Error probability 
(10~^) in realizing the photon transfer operation between 
A and _B as a function of the photon loss rate «: and CBJJ 
dephasing rate 
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For the pair of identical TLRs in FiglTfa), we introduce 
a single photon of frequency loq and it being in the leit or 
right TLR denotes the logic or 1 state for a single 
qubit. This is analogous to the conventional optical cav- 
ity mode representation of photonic qubit where the in- 
formation is encoded by which cavity the photon is in 0] . 
Notice for a dilution refrigerator temperature of 40mK, 
the thermal photon number in the TLRs is smaller than 
10~^° and thus the or 1 photon state for the TLRs is 
an excellent approximation. To effect arbitrary transfor- 
mations on this single qubit, we need to be able to shift 
the relative energies of the TLRs and transfer photons be- 
tween them, which implement the functionalities of phase 
shifters and beam splitters in optics. We realize this by 
coupling the TLRs capacitively to current biased Joseph- 
son junctions (CBJJ) as shown in Fig. [T](a). As the sim- 
plest Josephson qubit, CBJJ has the advantage that its 
level splitting can be easily adjusted by the bias current. 
Approximating the CBJJs as two-state systems with ad- 
justable energy splittings flc and [Sl, we can write 
the system Hamiltonian H = hji)Q[a)a + wb) + \hD.c<^l + 
^^ra^ + hgc[{a + b)a+ + (at + b^)a-] + hgr{b(j+ + 
[9|, where cr^'j^ are the Pauli matrices of the coupling 
and right CBJJ, a, b are the annihilation operators for 
photons in the two TLRs, and the coupling strengths 
gc,r = loqCctI \/2C{Cy+2C^, Cj'^ the capacitance 
of the couphng and right CBJJ. 

Since the CBJJ energies can be easily adjusted by tun- 
ing the bias current, we can control the interactions be- 
tween the TLRs and CBJJs. To transfer photons between 
the TLRs, we adjust the bias currents of the CBJJs to 
tune fir faraway from loq so the right CBJJ has no ef- 
fect. We further tune the coupling CBJJ close to reso- 
nance with luq and work in the dispersive region where 
the magnitude of detuning Ac — Qc — ujq is much greater 
than gc- Assuming the CBJJ was prepared in the ground 
state, its virtual excitation gives rise to the following ef- 
fective Hamiltonian for the TLRs 'i'l in the rotating frame 
defined by the uncoupled TLR Hamiltonian: 
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{ab^ +a)b) + ^{d)a + Pb). (1) 



Since there is only 1 photon in the system, a)a + w b = 1, 
energy shifts described in the second term in i?e// is a 
constant. The first exchange term implements the photon 
transfer operation. A photon can be transferred between 

2 

the two TLRs with a rate , which is about 20MHz 
for C} = 0.5pF, Cc = 23fF, and = 2GHz [l3,[ll|. 

To shift the relative energies of the TLRs, we tune 
the coupling CBJJ far off resonance and tune the right 
CBJJ into the dispersive region. Similarly, an effective 
Hamiltonian {hg^/Ar)b''b results which gives a relative 
phase when the photon is in the right TLR. 

We need to study the decoherence properties of our 
scheme to analyze its reliability. The photonic qubits 



have superb coherence and their life times are orders of 
magnitude longer than that of superconducting qubits. 
For TLRs fabricated on superconducting chips, a high 
quality factor of 10^ — lO'' has been demonstrated [12| . 
For TLR frequencies of tens of GHz, the photon loss rate 
k/27t can be as low as KHz. In contrast, the CBJJ has 
a short decoherence time, and we assume its dephasing 
rate r2/27r «lMHz. The CBJJ's decay rate from the 
excited state ri/27r is on the order of 0.1 MHz. 

A major advantage of our scheme is that the rela- 
tively lossy CBJJ does not damp the coherence of the 
photonic qubits much since it is only virtually excited. 
The CBJJ's decay from the virtually excited state in- 
creases the photon's loss rate by {gc,r/ Ac,r)'^^i, which 
is not a concern since {gc,r/ Ac,r)'^^i is no greater than 
K. To study the effect of the CBJJ's dephasing rate 
we model the dephasing effect as the result of a random 
fluctuation S„ in the CBJJ's energy splitting. This intro- 
duces an uncertainty in the detuning during for instance 
a photon transfer operation, Ac = f2c — ~* + ^n- 
Therefore, the system will have a random Hamiltonian 
Hnoise = —^{gd 5n{cL^b + ow) in addition to that in 
Eq. 11]). Assuming the distribution of 5„ is Gaussian, we 
can estimate the photonic qubit 's decoherence time due 
to Hnoise by using the free induction decay function [l^l ■ 
Since the free inductor decay function is determined by 
the spectral density of 5m which in turn is related to the 
CBJJ's dephasing rate it can be estimated that the 
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system's dephasing rate is no greater than 2(-^)r2. 

Following the quantum theory of damping, we now cal- 
culate the gate error of a photon transfer operation un- 
der the influence of cavity loss and CBJJ dephasing. us- 
ing the Master equation for the qubit 's density matrix p, 
dp/dt = —i[Heff,p] + K[apa^ — ^d^dp— ^pd^d] + K[bpb^ — 

fb^bp - ip6tS] + 2(f^)2r2[(dt6 + db^)p{d^b + aSt) - p]. 

The gate error probability of a single qubit bit flip is plot- 
ted in Fig. [Ijb) as a function of k and r2. The result 
indicates that, for already demonstrated r2/27r = IMHz 
and K/2n = lOkHz the gate error is on the order of 

The manipulations demonstrated so far perform linear 
optics. We still need a mechanism to induce interactions 
between photons. This is a major difficulty in conven- 
tional optics. However, at microwave frequencies, we can 
take advantage of the strong nonlinearities in Josephson 
devices to interact photons. 

We consider the low current biased 4-junction SQUID 
(FJS) device [l^ in Fig. [l{a). The two small identical 
SQUIDs are coupled inductively to TLR C of length / at 
positions ±//4. (This does not mean that the FJS must 
extend to a length of 1/2 because the TLR can be layed 
out in a zig-zag fashion.) Since / is much larger than 
the dimension of the FJS, we can adopt the long wave 
approximation and use the TLR current at the SQUIDs' 
locations in calculating the SQUIDs' flux bias. At the 
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FIG. 2: (a) The 4-junction SQUID to interact photons, (b) 
Dependence of the error probabihty (10"'') of a controlled 
phase gate on the ratio between photon transfer rate A and 
uncertainty in photon energy shift Us . (c) The setup for con- 
trolled photon interaction. The circuit in the dashed-line box 
is that in (a). The TLRs are coupled to CBJJs (not shown) 
so that single bit operations including photon transfer can be 
performed. 

two coupling points, the TLR currents are the largest in 
magnitude and opposite in direction. The main loop is 
coupled to TLR D at I /A. 

Assuming there is no other external flux bias, the small 
SQUIDs and big loop are biased by the TLR currents 
Ic — T'iIco{c — c') and Id = —ilrn){d — d^), where 
Icq = ^Jhwcl Lc and Ijjq = ^Jhwd/ Ld the zero point 
current fluctuations in the TLRs, Lc,d the inductance of 
the TLRs and c, d the annihilation operators for photons 
in C and D. We can work out the system's Hamiltonian, 
H = HTLR+Hpjs+H^nt, where Htlr = h{ujc+^s)c'^^+ 
h{ujd + Lo'^)dU, Hp. IS - - ^^-^-4^5cos(/., and 

Hint = hUmtC^ Cd'' d = hUinthcfld- (2) 

In these equations, ~ (2e)^/4(Cj + Cg) and Ej ~ 
hlc/2e are the charging energy and Josephson energy of 
the junctions, where Cj and /c are the junctions' ca- 
pacitance and critical current, (p is the average phase 
of the 4 junctions determined by the low bias current 
Jb « of the FJS. The frequency shift = -2E'j{(j>^xl + 
xlxDiK '^s = -2S5(0'x^ + xlxD/fi where Xc = 
wMcIco/{7rLsIc + $o), Xd = t:MdIdo/[tt{Ls + Ll)Ic + 
$o], Ls and Ll are the self-inductances of the small 
SQUIDs and the circuit loop. To simplify the expres- 
sions, we set X — Xc — Xd and denote the photon fre- 
quency shift ojs = ojg — ujg. The photon interaction 
strength ojmt — — iEjx'^ cos (j)/h. In deriving the system 
Hamiltonian, we have used the rotating wave approxima- 
tion and dropped terms that will be oscillating fast in the 



rotating frame deflned by Htlr- We have also dropped 
terms involving creation and annihilation of two photons. 
These terms have no effect since there is no more than 1 
photon in the TLRs in our scheme. 

We operate with a low bias current « for the FJS 
so that {coscf)) is large and the FJS' energy splitting is 
far away from the frequencies of the TLR photons. Thus, 
the FJS will not be excited by the TLR photons and they 
hardly get entangled. The FJS then acts as a "nonlinear 
medium" and Eq. ([2|) describes the interaction between 
photons in C and D modulated by the FJS' phase. For 
Ic = SO^iA, Ls ~ lOpH, and Mc « 80pH [11], the photon 
interaction strength ojint ~ IMHz, much greater than the 
photon loss rate. Unfortunately, there are difficulties in 
using this interaction for quantum computing. First, 4> 
has fluctuations in it due to the FJS' charging energy and 
thus the interaction strength is not a constant. Also, it 
is not easy to turn off the interaction. Tuning (p close to 
7r/2 requires biasing the FJS close to its critical current 
which makes the system unstable. The uncertainty in (p 
grows too. 

To have the photons interact only when needed, we 
use a setup shown in Fig. [2] (c). Here TLRs A, B and 
E, F are two qubits with photons being in A and F 
representing their logic state. When both qubits are 
in the 1 state, we can use the photon transfer operation 
discussed earlier to transfer the photons from B and E to 
the auxiliary TLRs C and D whose frequencies are made 
different than that of the qubit TLRs by uJs to account 
for their energy shifts. Once the photons are in C and D 
they can interact due to coupling to the FJS. Afterwards, 
we transfer them back to B and E. 

To stabilize the FJS' phase, we shunt its junctions 
with large capacitances Cg as shown in Fig. EJa). 
At low bias currents the FJS' behavior can be very 
well approximated by that of a harmonic oscillator 
and the distribution of (j) is given by its ground state 
wave function ^a/y^exp]— — 0o)^/2], where a — 
\/4:E^ cos (po /Ec and (po = {4>) = arcsiug^. If we 
choose a total capacitance Cj + Cg ~ 20pF, the rela- 
tive uncertainty 5{ujint)/oJint ~ 10^*. Such a small error 
is not a concern for the photon interaction term. How- 
ever uncertainties in the photon energy shift terms hujs 
can be comparable to hiOint and can cause large errors. 

We employ a two-phase technique in the spirit of spin- 
echo to address this problem. In phase 1, we first 
do a photon transfer operation between _B, C and i?, 
D with a speed relatively fast compared to ujint and 
SlOs = —2Ejx'^S{<p'^)/h, the uncertainty in the photon 
frequency shift. We then wait for a desired interaction 
time t = Tr/ujint after which we do another photon trans- 
fer between B, C and E, D. In phase 2, we first perform 
a bit flip for the 2 qubits, in other words do a photon 
transfer operation between A, B and E, F. We then re- 
peat phase 1. At the end, we perform a bit flip on the 



two qubits again. In this process, depending on their ini- 
tial states the qubits will acquire the same random phase 
due to 5ujs in either phase 1 or 2, thus removing the effect 
of the randomness in the photon energy shifts. The end 
result is a TT phase shift on the 2-qubit states if they are 
both in or 1 initially. This is equivalent to a controlled 
phase gate and it enables universal quantum computing 
in combination with the single qubit operations. 

If the photon transfer operation between _B, C and 
D was perfect, the controlled phase gate would be 
exact, phase shift in C and D could be eliminated com- 
pletely. However, since the photons in C and D will 
interact with the FJS even during the photon transfer, 
our control phase gate will have errors. This can be seen 
by examining the system Hamiltonian during the photon 
transfer (in the rotating frame) H — hXi,c{b''c -I- 6c^) -I- 
hXde {d'fe + de'')-2E°jX^(P'^{c^c + S d) - hio^^tc^ cS d. The 
first two terms are used for the photon transfer opera- 
tion, however the remaining terms cannot be turned off 
making the photon transfer imperfect. Obviously, the fi- 
delity of our controlled phase gate will be improved by 
making the photon transfer frequency Abe and Xde large 
compared to 5u)s and ujint- We numerically studied our 
control phase gate using the full Hamiltonian and plot- 
ted the gate error in Fig [2] (b). We set A = Abe = A^e- 
For our choice of system parameters, X/Scug ~ 20 and the 
gate error is on the order of 10""^. 

Our microelectronic system is easily scalable as shown 
in Fig. [5] (c). We can extend the setup for the control 
phase gate in both ends to integrate many TLR qubits on 
the same chip with an FJS between each pair of qubits. 
This is a Id architecture with controllable interactions 
between adjacent qubits that can be scaled to a large 
number of qubits. 

In order to perform photonic qubit quantum comput- 
ing, we still need to be able to generate and detect single 
photons. Photon generation on superconducting chip has 
been demonstrated experimentally 17|, |l8|, ll9| ■ For pho- 
ton detection 2^, 21 1 , we again consider a CB J J coupled 
to a TLR of frequency ujq as shown in Fig. [3] (a). The 
CBJJ is prepared in the ground state \g) in the well of its 
washboard potential. We also make use of an unstable 
excited state |e) where the CBJJ can tunnel to the volt- 
age state with a large rate F. By adjusting the CBJJ's 
bias current, we can tune the CBJJ in resonance with ujq. 
The CBJJ will then be excited by the TLR photon to |e). 
When it escapes from |e), a detectable voltage appears 
across the CBJJ. 

Though an easy and reliable method, our scheme may 
fail to detect a photon in the TLR due to the pho- 
ton decay and the CBJJ's intra- well decay and deco- 
herence. The TLR photon may decay before being de- 
tected by the CBJJ. The CBJJ's intra-well decay from 
|e) to \g) and its finite decoherence time are concerns 
too. To study the influence of the photon loss rate 
and CBJJ's intra-well decay and decoherence on the ef- 




FIG. 3: (a) The photon detection scheme, (b) The depen- 
dence of detector efEciency(l — 10~^) on the ratio between the 
escape rate F and photon loss rate k. The coupling strength 
gtd/2Tv = lOOMHz, photon loss rate K/2-n- = lOkHz, CBJJ 
decay and dephasing rate 7T/27r = lOOkHz, 7<^/27r — IMHz. 



ficiency of our photon detector, we model it as a 3- 
state system shown in Fig. [3] (a), where |/) repre- 
sents the voltage state. We use the Master equation 
dp/dt = —i[H,p\ -\- Cp. Here, p is the density matrix 
of the system, the system Hamiltonian H = Sa^a -\- 
9td{d'^<^ge + daeg), the detuning S ~ luq — p, p, the fre- 
quency difference between |e) and \g). The Liouvillian 
Cp = ^C[d]p+^C[a,f]p+^C[aeg]p+^l:^=g,eC[\^){^\]p, 
where £[6]p = 20^6+ - QtOp _ pOtQ. k is the de- 
cay rate of the photon in the TLR, 77^ and 7^ are the 
intra-well decay rate and dephasing rate of the CBJJ, 
= l«>0'l fo'' = 5>e,/, and cr^ = |e)(e| - \g){g\. 

Assuming initially there is a photon in the TLR and 
the CBJJ is in \g), we plot the detecting efficiency (the 
probability the CBJJ ends up in |/)) in Fig. [3] (b) as 
a function of F/k. As can be seen, the efhciency is 
high even for moderately large escaping rate F. For 
F/27r = 20MHz ^ and k/2tt = lOkHz, the detection 
efficiency is above 99%. Also, it is demonstrated in our 
simulation that the influence of to the detection ef- 
ficiency is minor and thus the detecting CBJJ docs not 
need to have long decoherence times. 

In summary, we have shown that, by using TLR mi- 
crowave photons as qubits and Josephson devices as opti- 
cal devices, a superconducting chip provides an ideal im- 
plementation for fully integrated photonic qubit quantum 
computing. Thanks to our careful design, high gate fi- 
delities can be achieved and thus our scheme is a realistic 
approach. Since our system is based on existing mature 
technologies, fast experimental progress can be expected 
to bring integrated photonic qubit quantum computing 
to reality. The novel idea of using on-chip microwave pho- 
tons as qubits also opens the possibility of investigating 
many interesting optical quantum effects in an integrated 
system. 
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